Proteins: Structure, Function, and Bioinformatics

Summary: Prot2Surf is a software tool designed for the characterization and prediction of protein association to surfaces. In this application note, Prot2Surf was tested using catalytic domains of the lytic polysaccharide monooxygenases (LPMOs), interacting with native surfaces. The results show that the software can efficiently analyze key binding features, including protein-surface distances, distances between catalytically reactive atoms, and the orientation angle between surface chains and the protein. These features are essential for distinguishing productive binding poses in these protein-surface systems and for understanding interaction patterns that provide guidance on protein engineering. Prot2Surf performs these analyses within seconds to a few minutes, providing a fast and accessible framework to post-process and characterize protein-surface encounter complexes. Availability and implementation: Prot2Surf, which is written in Fortran90, is documented and freely available as open source on GitHub: https://github.com/TUNNELING-GROUP/Prot2Surf. In order to run Prot2Surf, users should also install the SDA software package which is freely available at https://www.h-its.org/downloads/sda7/.

Quorum sensing in Gram-positive bacteria commonly relies on posttranslationally modified peptide pheromones. In Bacillus subtilis, the prenyltransferase ComQ catalyzes tryptophan prenylation of the quorum-sensing peptide ComX, but the structural basis of this unique peptide modification has remained unclear. Here we identified a previously uncharacterized ComQ homolog, StheQ, and its cognate peptide substrate, StheX, from Sphaerobacter thermophilus and investigated their structural and functional relationship. Liquid chromatography-tandem mass spectrometry (LC-MS/MS) analysis demonstrated that StheQ catalyzes prenylation of the tryptophan residue located second from the C-terminus of StheX. Crystal structures of apo StheQ and its complexes with a farnesyl pyrophosphate analog revealed that StheQ adopts the all--helical fold of the trans-isoprenyl diphosphate synthase (IPPS) superfamily while possessing an active-site architecture adapted for peptide-based indole prenylation. The structures identified a single Mg2+-binding site associated with the first aspartic acid-rich motif and showed no evidence for metal coordination at the pseudo-second aspartic acid-rich motif. Site-directed mutagenesis, complex formation assays, and docking analyses identified a peptide-binding pocket adjacent to the active site and suggested that N215 contributes to productive positioning of the acceptor tryptophan. These findings establish the structural basis for peptide prenylation by a ComQ-family enzyme, providing insight into the evolution of peptide-based indole prenylation within the IPPS superfamily, and support the view that ComQ-family enzymes constitute a distinct functional branch specialized for peptide modification.

The rotary motor F1-ATPase has been extensively studied as a model molecular machine, yet rational engineering of its catalytic activity remains challenging because ATP hydrolysis is regulated by long-range intersubunit allostery and large conformational transitions. Here, we developed a homolog-guided engineering strategy to increase the maximum rotation rate of the thermophilic Bacillus PS3 F1-ATPase (TF1). Candidate mutation sites were first identified by comparing TF1 with the homologous enzymes bovine mitochondrial F1 (bMF1) and Paracoccus denitrificans F1 (PdF1), both of which exhibit higher maximum rotation rates than TF1. Systematic exploration of these sites identified four activity-enhancing hotspots, followed by focused hotspot exploration and machine-learning-assisted prioritization of combinatorial mutants. The best mutant, TF1({beta}Y313L/{beta}E332S), exhibited a 1.8-fold higher maximum rotation rate than TF1(WT) while retaining its functional thermostability. Interestingly, activity-enhancing substitutions were not limited to the residues conserved in both bMF1 and PdF1, indicating that the bMF1-PdF1 consensus substitutions effectively identify activity-enhancing hotspots rather than uniquely defining the optimal amino acid. Machine-learning-assisted exploration efficiently prioritized highly active mutants, although the predictive performance was limited by the relatively small training dataset and epistatic interactions among mutations. Kinetic and structural comparisons further provided mechanistic insights into the enhanced catalytic activity of the engineered mutant. Together, these results establish a practical strategy for engineering complex molecular motors by combining homolog-guided hotspot identification with focused hotspot exploration.

The folding thermodynamics of proteins are dominated by two opposing forces, the loss in backbone entropy and the packing of hydrophobic groups. The same forces are major contributors to the extension thermodynamics of elastic proteins with the distinction that both processes act in concert, favoring the higher chain and solvent entropy of a relaxed conformation. The relative entropic contributions specify the recoil mechanism; human elastin recoil is primarily driven by hydrophobic forces, whereas fly resilin has a rubber-like mechanism driven by backbone entropy. Despite the importance of elastic proteins to tissue biomechanics, few have been identified, let alone characterized to the same extent as elastin and resilin. We develop a thermodynamic framework that maps proteins by sequence-derived estimates of extension-induced backbone and solvent entropy changes. Putative elastic proteins are proposed and classified by recoil mechanism based on estimated thermodynamic features. Proteins that map to elastic regions are overrepresented by the skin proteome. The set of predicted elastic domains is further extended by incorporating sequence context embedded in protein language models. Protein domains with distinct thermodynamic recoil mechanisms cluster on the latent space manifold. Some of these domains are anticipated to have roles within molecular machines, expanding the scope of elastic protein function beyond mechanical materials like elastin and resilin.

Human cytochrome P450 2C9 (CYP2C9) is a hepatic microsomal enzyme involved in the oxidative metabolism of clinically important drugs, but the structural organization of its oligomeric assemblies outside crystallographic packing environments remains poorly understood. Here, we report the cryo-EM structure of human CYP2C9 determined under aqueous, membrane-free conditions at 3.31 Angstrom resolution. The structure reveals a C2-symmetric hexameric assembly organized as a dimer of trimers. Individual protomers retain the conserved P450 fold and heme-binding architecture observed in previously reported crystal structures, indicating that assembly formation does not substantially perturb the catalytic core. The hexamer is stabilized by defined intra-trimer interfaces involving the N-terminal region and residues around Trp212 and Phe482, together with inter-trimer interfaces involving Leu71 and the 220-227 loop. These interfaces are distinct from the crystal packing contacts observed in CYP2C9 crystal structures, demonstrating that the assembly is not a simple recapitulation of crystallographic packing. Notably, the inter-trimer interface is located near the FG-loop-containing surface previously implicated in membrane association. This suggests that the observed hexamer may represent a membrane-free association of two trimers through membrane-related surfaces, whereas the trimeric arrangement itself may be compatible with membrane-associated organization. The structure therefore provides a framework for investigating how trimer formation, membrane interaction and local conformational changes in the FG-loop region may influence CYP2C9 function.

BACH1 is a heme-regulated basic-leucine-zipper containing transcriptional repressor that binds its DNA recognition elements as a heterodimer with MAFK. Heme-binding is thought to be mediated by several Cys-Proline (CP) motifs and this results in dissociation of the heterodimer from DNA. The mechanism of heme-binding and heme-mediated DNA dissociation remains unresolved. We have used UV-visible spectroscopy, 2D-NMR and DNA-binding assays to explore both heme-binding and DNA dissociation of a minimal BACH1 construct containing 2 CP motifs (C492(CP5) and C646(CP6)) flanking the DNA-binding domain. We find that heme is able to bind to both CP motifs, but also to other non-CP cysteines and histidines in the construct. Using NMR spectroscopy, we identify a structured binding pocket in which heme interacts with both C646(CP6) and Cys621. However, DNA-binding assays show that C646(CP6) is not required for heme-mediated DNA dissociation of the BACH1:MAFK heterodimer. Using UV-visible spectroscopy we show that C492(CP5) also recruits heme with a second ligand, a conserved histidine, His559, in the BACH1 DNA-recognition helix. Mutation of C492(CP5) reduces but does not abolish heme-mediated dissociation from DNA. Our findings suggest a mechanism for heme-binding to BACH1 and heme-mediated dissociation from DNA.

Aptamers are structured nucleic acid ligands capable of high affinity, high specificity molecular recognition generated using variations of the SELEX (Systematic Evolution of Ligands by Exponential Enrichment) process. However, SELEX often produces sequences that enrich yet may lack binding efficacy. We propose a measure called the Ruggedness Composite Index (RCI) along with a method for computing it, that can be used to distinguish binding-competent ('active') aptamers from weak or non-binding ('inactive') aptamers. Given a set of aptamers, RCI incorporates information on their fragmentation (landscape partitioning), basin entropy (metastable state distribution), cumulative density irregularity (non-uniform occupancy), and structural energy correlation length (structure-energy coupling scale). We test whether secondary-structure folding energy landscape topology distinguishes active from inactive aptamers using a multiscale level set framework across six datasets. Active aptamers show lower RCI values and occupy smoother, funnel-like conformational spaces, while inactive aptamers show higher RCI values, reflecting fragmented, high-entropy landscapes. By contrast, classical thermodynamic features, such as minimum free energy, show limited discrimination between active and inactive aptamers. In all datasets, sequences that exhibit enrichment which is not monotonic but lack specificity exhibit elevated ruggedness, indicating landscape topology can predict non-specific enrichment. These results indicate that folding landscape organization can be used as a predictor of aptamer activity and establish RCI as a simple, mechanistically interpretable measure for improving candidate prioritization, especially in therapeutic aptamer discovery.

Alzheimer's disease is characterized by the accumulation and aggregation of amyloid-{beta}(A{beta}), but the molecular mechanisms linking environmental and infectious factors to A$\beta$ conformational changes remain incompletely understood. Herpes simplex virus type 1 (HSV-1) has been proposed as a potential contributor to AD pathology, and interactions between the viral glycoprotein B (gB) and A$\beta$ may influence the conformational behaviour of the peptide. Molecular dynamics (MD) simulations provide atomic-scale information on such interactions, but conventional structural descriptors may not fully capture changes in the organization of residue interaction networks. Here, we introduce a graph-geometric framework based on Forman-Ricci curvature to characterize the evolution of residue interaction networks during MD simulations. Each simulation frame is represented as a residue interaction graph based on C--C contacts, and residue-wise curvature profiles are analysed across time. We apply the framework to A{beta}1-42 in isolation and in complex with HSV-1 gB. Conventional MD analyses indicate stable association of the simulated complex, favourable interaction energetics, and conformational changes in A{beta}, including a transition from -helical structure toward {beta}-turn-rich conformations over the simulated timescale. Forman-Ricci curvature reveals pronounced and spatially localized remodelling of the A{beta} residue interaction network in the complex, with the strongest changes concentrated in the C-terminal region. These regions also exhibit reduced temporal curvature fluctuations and progressively distinct geometric behaviour throughout the simulation. Hierarchical clustering further identifies cooperative groups of residues with coordinated curvature dynamics, including a prominent C-terminal domain. Together, these results demonstrate that Forman-Ricci curvature provides a complementary description of biomolecular dynamics by capturing changes in the geometric organization of residue interaction networks that are not directly represented by conventional structural descriptors. The framework provides a general computational approach for studying network-level structural remodelling in protein molecular dynamics and offers a quantitative perspective on the conformational consequences of HSV-1 gB--A{beta} association.

Amyloid fibrils are highly ordered protein aggregates characterized by a conserved cross-{beta}-sheet architecture despite originating from structurally diverse precursor proteins. Growing evidence suggests that interactions between different amyloidogenic proteins can modulate aggregation pathways through heterologous cross-seeding; however, the influence of seed polymorphism on the structure and biological properties of cross-seeded fibrils remains poorly understood. Here, we investigated the cross-seeding of native human insulin by two structurally distinct polymorphs of hen egg-white lysozyme (HEWL): flexible fibrils (FFs) and rigid fibrils (RFs). Native insulin remained stable under physiological conditions and underwent spontaneous fibrillation only under acidic conditions. In contrast, both HEWL polymorphs efficiently induced insulin aggregation at physiological pH, bypassing the nucleation barrier. Thioflavin T fluorescence, circular dichroism spectroscopy, and transmission electron microscopy revealed that lysozyme FFs templated the formation of insulin flexible fibrils (IFFs), whereas lysozyme RFs produced insulin rigid fibrils (IRFs), demonstrating that the structural characteristics of the parental HEWL polymorphs were propagated during heterologous cross-seeding. The toxicity of the resulting insulin fibrils was evaluated in SH-SY5Y neuronal cells and CCF-STTG1 astrocytes. IFFs exhibited minimal cytotoxicity and only subtle morphological alterations, whereas IRFs caused modest reductions in cell viability accompanied by more pronounced cellular damage. These findings demonstrate that the structural polymorphism of HEWL fibrils governs both the architecture and biological activity of cross-seeded insulin fibrils, highlighting amyloid polymorphism as an important determinant of heterologous amyloid propagation and a potential design principle for engineering functional amyloid-based biomaterials and protein delivery platforms.

Infectious bronchitis virus (IBV) is a member of the Gammacoronavirus genus responsible for respiratory illness and weakened eggshells in infected chickens, adversely impacting the poultry industry. Escaping innate immune detection during infection is crucial for coronavirus proliferation in the host. The production of double-stranded RNA during coronavirus replication triggers innate immune sensors to create an antiviral state within infected cells. To counter this response, coronaviruses employ nonstructural protein 15 (nsp15) endoribonuclease to degrade double-stranded RNA. Here, we use cryo-electron microscopy and biochemistry to characterize IBV nsp15 interactions with RNA. While the overall structure and active site of IBV nsp15 strongly resemble previous studies of nsp15 from other coronaviral genera, we note that double-stranded RNA contacts several non-conserved residues peripheral to the enzyme active site. Our data show that these residue positions can have strong impacts on RNA cleavage suggesting unique solutions for RNA engagement across coronavirus species. We also demonstrate a preference for IBV nsp15 to cleave double-stranded RNA over single-stranded RNA and observe nsp15 hexamers with two double-stranded RNAs bound simultaneously. This study reinforces the need to study diverse coronavirus species to identify distinct viral enzyme characteristics.

Biomolecular condensates formed by intrinsically disordered proteins require molecular models that accurately describe proteins in both dilute solution and condensed phases. Explicit-solvent coarse-grained models offer an attractive balance between chemical resolution and computational efficiency. Yet, it remains unclear whether improving dilute-state properties is sufficient to obtain an accurate description of condensates. Here, we address this question by introducing minimal modifications to the Martini 3 force field that combine recent advances in bonded interactions with refined protein-water interactions and strengthened glycine self-interactions, while preserving the underlying chemical transferability of the model. The resulting model substantially improves the description of single-chain conformations across a diverse benchmark of disordered proteins. We then investigate phase separation of the well-characterized low-complexity domain of heterogeneous nuclear ribonucleoprotein A1 and its sequence variants. The model reproduces several key physicochemical properties of biomolecular condensates, including chain expansion in the dense phase, sequence-dependent intermolecular contacts, protein diffusion and its relation to single-chain dimensions, and hydration, while revealing quantitative limitations in condensate density, phase equilibria, and ion partitioning. Our results show that improving dilute-state behaviour translates into a better description of condensed-phase properties, including condensate density, but is not sufficient to quantitatively reproduce the equilibrium between the dilute and dense phases.

Scoring biomolecular complexes is central to structure assessment and drug discovery, yet the complexes themselves vary widely in pose, size, and molecular composition. A scoring function tuned for one interaction type rarely carries over to another, and most existing methods compound the problem by leaning heavily on task-specific labels. We introduce OmniScore, a universal structure-based framework that learns a shared geometry-aware representation of complexes once and then adapts it to downstream scoring through lightweight task-specific heads. OmniScore couples a graph view and a sequence view of each structure, encodes its three-dimensional geometry, and compresses representations into a compact latent space that a reconstruction module and prediction heads can reuse. We pretrain this backbone on diverse datasets including complexes, monomers, and small molecules with complementary objectives: coordinate recovery, correcting corrupted input tokens, predicting molecular identity, and grounding the representation in structure-level physical quantities. Across the evaluated benchmarks, OmniScore gave the best antibody-antigen and nanobody-antigen quality assessment on all reported metrics compared to state-of-the-art baselines. Its frozen residue embeddings matched the state-of-the-art protein-tokenization method with an average functional-site accuracy of 71.8% on a standard residue-level benchmark. On protein-ligand scoring and ranking benchmarks, it performed on par with methods built specifically for that single task. These results suggest that geometry-aware pretraining can provide a reusable scoring backbone for tasks that depend on interfacial and residue-level structure, within the evaluated settings.

Deep learning-based protein structure prediction methods that leverage evolutionary information from multiple sequence alignments (MSAs), exemplified by AlphaFold2, have achieved remarkable accuracy. However, existing methods still struggle to predict challenging proteins, particularly those with novel folds or limited evolutionary information, and to recover alternative conformational states. Here we show that structure prediction models trained under different MSA-depth distributions corresponding to different levels of evolutionary information exhibit complementary generalization behaviors, and that a model trained on a mixture of these distributions can combine their complementary generalization strengths. Building on this insight, we developed ProtMonomer, a deep learning framework trained on MSA-depth distributions representing a broad range of evolutionary information levels to improve structure prediction. Across benchmarks comprising CASP15 targets, non-redundant experimentally determined structures, orphan proteins, and short peptides, ProtMonomer performed comparably to or better than leading methods, including AlphaFold2 and AlphaFold3, with particularly strong performance on challenging targets. For fold-switching proteins, ProtMonomer also recovered alternative conformational states more accurately than AlphaFold2 and AlphaFold3 across diverse homologous sequence sampling strategies. In addition to improving predictive accuracy, ProtMonomer substantially reduced inference cost through an efficient architecture, enabling high-throughput applications. Together, these findings provide insights into the generalization of evolution-informed structure prediction models and support ProtMonomer as an accurate and efficient framework for protein structure prediction.

Linker chemistry and conformation are central determinants of PROTAC activity, shaping ternary-complex geometry, cooperativity, target-lysine presentation and cellular permeability. Existing linker generators often lack explicit control over linker flexibility, require predefined attachment sites and linker lengths, or produce structures that demand substantial geometric correction, limiting their utility in practical PROTAC design. Here we introduce FlexiTAC, a Bayesian flow network that jointly generates linker atom types and coordinates from the warhead and E3-ligase-ligand contexts. We also assemble PROTAC-3D, a quality-controlled collection of 63,554 component-resolved PROTAC structures for model training, and PROTAC-Bench, which covers molecular quality, fragment preservation, geometric fidelity, conformational stability, fragment awareness, rediscovery and sampling efficiency. Compared to the best 3D baseline models, FlexiTAC improves validity by 12.0-12.7% and achieves the highest PoseBusters pass rate of 79.5%-80.0%. A differentiable guidance module shifted generated linkers along a conformational ensemble-derived rigidity axis without retraining the generator. In silico case studies further show that the model can accept crystal-derived, redocked or predicted structural inputs. Together, FlexiTAC, PROTAC-3D and PROTAC-Bench establish an integrated and reproducible framework for data-driven PROTAC linker design, combining controllable structure-conditioned generation with standardized training data and evaluation protocols. This framework expands the linker chemical and conformational space accessible to computational exploration, provides a foundation for future method development and enables the systematic generation of structure-conditioned linker designs with tunable conformational flexibility.

The binding affinities of a chimeric analog of a riboflavin derivative linked to biotin, (6- (7,8-dimethyl-2,4-dioxo-3,4-dihydrobenzo[g]pteridin-10(2H)-yl)hexyl 5-((3aS,4S,6aR)-2- oxohexahydro-1H-thieno[3,4-d]imidazol-4-yl)pentanoate), referred to as C6-Rf-biotin-tag, to the riboflavin binding retain or streptavidin are in the M range, 1.29 {+/-} 0.277 and 3.00 {+/-} 0.459, respectively. These values suggest that C6-Rf-biotin-tag has potential applications in diagnostic assay and labelling target flavin binding proteins. The C6-Rf-biotin-tag which was characterized with respect to physical and biochemical properties retains UV/Vis spectroscopic and fluorescence behavior similar to riboflavin.

Bacteriophages can only be understood through their interactions with bacterial hosts. As environmental sequencing efforts expanded, the number of available phage genome sequences has exploded, yet the vast majority of these sequences lack host information. Predicting the host of a newly observed phage is therefore a key challenge in virology. Several computational tools can predict phage-host relationships from genomic data, but they share notable limitations: (1) the number of different hosts that can be predicted remains relatively restricted; (2) tools tend to assign confident host predictions to non-viral input sequences; and (3) most tools have a trade-off between accuracy and speed. Here we present PhageTransformer (PT), a deep learning model for phage-host prediction that addresses these limitations. We benchmark PT against existing tools on 3,881 independent phage-host pairs from GenBank and public HiC data, and demonstrate that it achieves competitive or superior prediction accuracy at greatly reduced runtime.

Obesity is a growing public health challenge, with body mass index (BMI) influenced by both genetic and environmental factors. While the role of direct genetic transmission is well established, evidence for genetic nurture effects, in which parental genotypes impact offspring through the environment, has remained mixed. This study investigates direct genetic transmission and genetic nurture effects on BMI across ages, using parent-offspring trios and pairs from the Dutch Lifelines cohort study (N = 18,897 offspring, aged 8 to 67 years). We leveraged the latest multi-ancestry BMI polygenic score (PGS) to construct transmitted (PGS-T) and non-transmitted (PGS-NT) polygenic scores, where PGS-NT consists of parental alleles not passed on to offspring and serves as a proxy for genetic nurture. Linear mixed models showed a large effect of PGS-T on offspring BMI (Beta = 0.416, p < 0.001), corresponding to a 1.85 kg/m2 increase per SD increase in PGS-T. PGS-NT had a small but significant effect (Beta = 0.026, p = 0.013), consistent with a genetic nurture effect accounting for approximately 6.6% of the effect of direct transmission. Parent-of-origin analyses showed that maternal PGS-NT effects were larger than paternal effects. PGS-T interactions with age indicated that direct transmission effects increased in childhood and stabilized in adulthood, while PGS-NT effects remained stable across age. Our findings suggest that direct genetic transmission is the dominant influence on BMI, while results are consistent with small genetic nurture effects that are driven by the maternal side.

Background: In Denmark, perioperative high-dose glucocorticoid treatment were step-wisely implemented for total hip arthroplasty (THA), total knee arthroplasty (TKA), and unicompartmental knee arthroplasty (UKA). We aimed to estimate the effect of a single high dose of glucocorticoids on opioid consumption following primary THA, TKA, and UKA. Methods: This was a prespecified analysis of a multicenter natural experiment using electronic health record data. We included elective THA, TKA, or UKA surgeries performed in Eastern Denmark from 2017-2025. At each center, surgeries before implementation of high-dose glucocorticoids served as controls, whereas surgeries after implementation comprised the intervention group. The primary outcome was the between-group difference in cumulative 0-24h opioid consumption, which included preemptive end-of-surgery doses. The predefined minimal important difference was set at 5 mg IV morphine equivalents. Secondary outcomes were maximum 0-10 numerical rating scale (NRS) pain score and incidence of opioid-related adverse events within 24 hours, hospital length of stay, and days alive and out of hospital at 30 days. Results: A total of 47,317 surgeries performed at nine centers were analyzed: 13,010 controls and 34,307 in the intervention group. During the study period, five centers implemented high-dose glucocorticoids for THA patients, two for TKA/UKA patients. High-dose glucocorticoids were administered to 6% of patients before implementation versus 92% after. High-dose glucocorticoids resulted in a mean reduction of 3.8 mg intravenous (IV) morphine equivalents (95% CI 3.3;4.3). The intervention also reduced the maximum 0-24h NRS pain score by 0.8 points (99% CI 0.7;0.9), but there was no difference in adverse events, length of stay, or days alive and out of hospital. Conclusions: Implementation of high-dose glucocorticoids reduced 0-24-hour opioid consumption by 3.8 mg IV morphine equivalents after elective hip and knee arthroplasty. This difference was below the prespecified minimal important difference threshold. Online registration: https://doi.org/10.1101/2025.11.11.25339982

Micronutrient deficiencies are prevalent in Madagascar, where diets rely heavily on starchy staples and access to animal-source foods is limited. Small dried fish (SDF) are widely available, yet their nutritional value and health risks remain poorly documented. We combined market surveys, taxonomic identification, and micronutrient and heavy metal analyses of nine SDF types collected along National Road 7. The samples encompassed 33 fish families, were dominated by small pelagic species (Clupeidae and Engraulidae), and were appreciated by consumers. A daily portion (5 g for infants; 10 g for young children and women of childbearing age) contributed substantially to Recommended Nutrient Intakes (RNIs). Across samples and groups, SDF were rich (>30% of RNI) in selenium and, for infants and young children, in calcium. All samples were a source of (>15% of RNI), or rich in, phosphorus, whereas iron contributions were more variable but often substantial. Several samples exceeded 100% of RNIs for selenium, calcium, iron, or manganese in infants and young children, and some were also sources of magnesium and, less frequently, zinc. Vitamin A was absent from sun-dried samples but detected in a smoked freshwater type. Heavy metal concentrations varied markedly, and portions of several types led to estimated exposures to inorganic arsenic or cadmium exceeding reference values, whereas freshwater species and some pelagic types showed a more favorable nutrition-risk balance. Overall, SDF are affordable, nutrient-dense foods with strong potential to alleviate micronutrient deficiencies in Madagascar, while highlighting the need for type-specific guidance to balance nutritional benefits and contamination risks.

Abstract Background: Vaginal steaming has gained popularity among women for reasons known best to them. However, the effects of vaginal steaming on vaginal lactobacilli levels remain poorly understood and understudied. This study investigated the prevalence, assessed differences in the presence of bacterial vaginosis (BV) among women who practiced vaginal steaming and those who do not and determined the factors associated with vaginal steaming among women attending Mulago Sexually Transmitted Infections (STI) clinic in Uganda. Methods: This study utilized a cross-sectional design to enroll 181 women aged 18 to 49 years who were systematically sampled at Mulago STI clinic. Interviews were conducted to obtain the demographic characteristics of the participants. Vaginal swabbing and Gram staining were done to attain lactobacilli counts by microscopy which were categorized using the Nugent score. Data were analyzed using Stata. The potential confounding effects of other variables on the relation between vaginal steaming and the presence of bacterial vaginosis as well as factors associated with vaginal steaming were assessed using modified Poisson regression. Results: Prevalence of vaginal steaming was 40.3%, (95% confidence interval (CI) 33.0% - 48.0%). There were 41.1% women who practiced vaginal steaming occasionally, 53.4% who used hot water having herbs and 78.1% who practiced vaginal steaming for medical reasons. There was no difference in the presence of bacterial vaginosis when women who practiced vaginal steaming were compared to those who did not (p-value = 0.286). Factors that were significantly associated with vaginal steaming included having experienced vaginal issues (aPR = 0.07, 95% CI 0.01 - 0.12, p value = < 0.001) and contraceptive use (aPR= 0.52, 95% CI 0.37 - 0.72, p value = 0.001). Conclusions: About 2 in every 5 women at Mulago STI clinic reported to have indulged in vaginal steaming. There was no difference in the presence of bacterial vaginosis when women who practiced vaginal steaming were compared to those who did not. Having experienced vaginal issues and contraceptive use were significantly associated with vaginal steaming among women at Mulago STI clinic, Uganda. The Ministry of Health of Uganda should establish targeted screening and treatment for bacterial vaginosis alongside other sexually transmitted infections.